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Abstract 

We calculate the ground state current densities for 2+1 dimensional 
free fermion theories with local, translationally invariant boundary states. 
Deformations of the bulk wave functions close to the edge and boundary 
states both may cause edge current divergencies, which have to cancel in 
realistic systems. This yields restrictions on the parameters of quantum 
field theories which can arise as low energy limits of solid state systems. 
Some degree of Lorentz invariance for boosts parallel to the boundary 
can be recovered, when the cutoff is removed. 



Introduction 



The prediction of a relativistically invariant quantum field theory in 2+1 
dimensions as low energy limit of a solid state system [1] was confirmed 
by investigation of graphene [3]. Modification of the graphene structure 
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may allow to obtain more general quantum field theories [2]. Since solid 
state systems typically have edges, the study of quantum field theories with 
boundaries acquired additional interest. An important tool in this con- 
text is the bulk-boundary OPE. In particular, some of its singularities yield 
divergent observables and constitute obstructions to any experimental re- 
alization. The coefficient function of the identity in such an OPE is just 
the corresponding vacuum expectation value in the presence of the bound- 
ary. Integrated expectation values of current operators can be interpreted 
as measurable charge or spin currents and should remain finite in realis- 
tic systems. Here we study this obstruction in the simplest case of free 
fermions, where vacuum expectation values can be interpreted as integrated 
contributions of the particles in the Dirac sea. 

It is of course well known that such effects can prevent the lifting of a 
quantum mechanical system to a quantum field theory, even without the 
presence of a boundary. In our case the Dirac sea of a single fermion yields 
a half- integral value of the Hall conductivity [4] , [5] , [6] . Consequently one 
needs an even number of fermions, though some of them may be shifted to 
infinite mass, which results in a Chern-Simons term instead. For bound- 
ary states we shall find more intricate obstructions, which depend on the 
numerical parameters of the boundary state. 



2 Boundary conditions from self-adjoint extensions 

We consider spinor fields tp acted upon by the Dirac operator 

d d 

(1) H := -iai- 102^ + ma^i 

ox oy 

on the open half plane {x = S M^j x > 0}, as in [7]. This operator 

is a closable symmetric operator on the domain T>{H) = Cg°(M+ x M, C^) 
of smooth functions with compact support vanishing in a neighbourhood 
of X = 0. The possible boundary conditions are classified by its self-adjoint 
extensions. The latter can be obtained by the standard von Neumann theory. 
Here we give a modification, which is less general but better adapted to the 
present problem. Since 'D{H) is dense in the Hilbert space T-L := L'^{M.^ x 
M) (8) C^, the image of an arbitrary vector ^p £ Ti in the domain of a self- 
adjoint extension is uniquely determined by the unextended operator, so the 
problem is reduced to the determination of such domains. Let H* be the 
adjoint of the unextended H. By von Neumann theory, every self-adjoint 
extension of H corresponds 1-1 to an isometry V : A^-f — t- between the 
defect spaces 

N± = {'il;£ V{H*)\ H*tP = ±i^np}, 
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where /x is a fixed but arbitrary positive number. Indeed, V gives rise to the 
self- adjoint extension Hy with domain 

(2) V{Hv)=V{H) + {I-V)N+. 

The choice ^ = 1 is standard, but we shah see that the hmit of large fi 
allows a more direct physical interpretation. Applied to (1), N± is given by 
k integrals over elements of the one-dimensional spaces 



where 



s 



±^i{k + Xl^ and X = ^/ ^i^ + + m^. 



For ^ — 7- oo the defect spaces consist of wave functions which become sup- 
ported on an arbitrarily small neighborhood of the boundary and factorize 
in a universal term exp(— ^x) and a spinor of the form 

9{y)[l^ with 9GL2(M). 

Thus isometries — t- become tensor products of a unitary map V : 
L2(M) ^ l2(R) with the map between one dimensional subspaces of C 
which takes (1,-1) to (1,1). From (2) we read off that any tp G T){Hy) 
satisfies 

V'iU=o = (1 - y)9, 

i^2\x=o = {l + V)g. 
If 1 is not an eigenvalue of V , then 

(3) 'i/'2U=o = «r V'iU=o, 

where F is the self-adjoint operator 

T = -i{l + V){l-V)-^. 

The physical necessity of such a boundary condition can be understood as 
follows. Consider the conserved current densities 

f{-K) = lim VXx')"^<^''V'(x),. 

One needs a boundary condition which implies J j^{0,y) dy = 0, and this 
indeed follows from (3). 

We only study local boundary conditions, for which j^(0, y) vanishes iden- 
tically and also impose translational invariance. Thus F has to commute 
with d/dy and with multiplication by functions of y. This implies that it is 
a real constant 7, and we obtain the boundary condition 

■02(0,2/) = i7 1/^1(0,1/). 
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which we shah use in the sequel. To include the case V = 1, we consider 7 
as an element of the projective real line, including the value 7 = 00. 

For any 7 the Hamiltonian commutes with the PT transformation ip[x,y) 1— t- 
a3ip*{x,—y). The space reflection {x,y) 1— )• {x,—y) lifts by ip{x,y) 1— )■ 
aiip{x, —y) to a map between systems with parameters (m, 7) and (— m, —7"^) 
Because of this duality we only need to consider positive m. Arbitrary val- 
ues of m will only be considered at the very end of the calculation. It turns 
out that the results allow a continuous extension to m = 0, which is the case 
relevant for graphene. 

The CPT map ip{x,y) 1— )■ a2ip{x, —y) relates the positive energy states 
for (m,7) to the negative energy states for (m,7~^). A different type of 
duality is given hy xp{x,y) ^ a24'{—x,y), (m,7) 1— )• (—m, 7"^), which relates 
boundary states of complementary half planes. 

3 Eigenfunctions 

The Dirac operator commutes with —id/dy and consequently with 



We denote the real eigenvalues of these hermitean operators by k and t^. 
We take I G M"*" for positive P and / = i\ with A S M"*" for negative P. The 
corresponding eigenvalues E of H satisfy E'^ = k'^ + P + ir? . 

According to the sign of the spectrum of H has two parts, which will 
be refered to as bulk and edge. Bulk eigenfunctions u/fc with given /c, 1? have 
the form 



with constant spinors x±- The equation -ffn^fc = Euik and the boundary 
condition reduce to {m — E, ±l — ik)x± = and (x+ + X-)2 = ^7(x+ + X-)i- 
For each / > there is a unique solution up to normalization, such that for 
the bulk part we just recover the spectrum of the Dirac operator on the 
whole plane, with a gap A = (— m,m). We normalize these eigenfunctions 
with respect to the measure dk dl j (27r^ ) and choose phases such that 



{idjdxf = H'^ + {d/dyf 



— m 



2 



uik{x, y) = exp(i/x + iky)x+ + exp(-i/x + iky)x 




E — m 



4E 
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where 

with 

k + il 



(5) Pik ■-- 



m — E 



Edge eigenfunctions [7] have the form Uk = exp(— Ax + iky)x- We nor- 
mahze them with respect to the measure dk/ir and choose phases such that 



The equation HUk = EUk yields 



and 



27 , 1-72 
E = - — '-^k + 



7^ + 1 7^ + 1 



Thus E, A are one- valued functions of k, and for any k one has a unique edge 
state whenever A > and no edge state otherwise. The latter inequality 
restricts the possible values of k, E. The limiting values for A = lie on the 
mass shell. When we put 

1-7 
V = sgn— — , 
1 + 7 

then positive energy half of the mass shell is touch for mr] = 1 and the 
negative energy half for mrj = —1. One finds energies within the spectral 
gap A of the bulk iff > 0. 

The linear dispersion relation for the boundary equation implies that edge 
states travel at the fixed velocity 



27 
1 + 7 

Note that v^dge determines 7 up to the replacement of 7 by 7^"*^, which is 
the one induced by CPT. The latter transformation changes the signature 
Tj. Altogether, the boundary condition can be characterized by Vedge and ry, 
except for the fact that rj is undefined in the limiting cases 7 = ±1. One 
observes that 
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1) for vanishing edge velocity, i.e. for 7 = and 7 = 00, edge states have 
E = m and E = —m, respectively, independently of k. This case is 
probably unphysical, since all edge states have the same energy and 
the thermodynamic partition function diverges for a system of infinite 
width. Since the width of the edge excitation increases linearly with 
k, the thermodynamic properties of a system of a finite system would 
depend on the sample size. 

2) The edge velocity is equal to the maximal bulk velocity for the CPT 
invariant boundary conditions 7 = ±1. These yield E = jk and 
A = 7m independently of k. For 7m < there are no edge states, for 
7m > all edge states have the same dependence on x. This yields 
particular divergences, so we will tacitly exclude 7 = ±1 in most of 
the subsequent calculations. 

In (k, ii^)-space the half-line of edge states is tangent to the half-hyperboloid 
given by Z > and E"^ = k'^ + m"^ + P at I = X = 0. For 7^ < 1 this happens 
at positive energy and for 7^ > 1 at negative energy. The tangency has 
consequences for the behaviour at large distance from the boundary, as we 
shall see later. 

The portion of edge energies within A gives rise to an edge current parallel 

2 

to the boundary, with edge-conductivity equal to (in units of ^) [7] 



This applies whenever some subinterval of the gap is actually occupied, 
independently of the size of this subinterval. 

4 The current densities 

When no boundary exists, Fermi energies within the gap A obviously yield 
the same theory. This means that any pair (7, Ep) with Ep G A yields 
a boundary state of the standard Lorentz invariant free fermion theory. 
We want to calculate the corresponding expectation values of the current 
densities. The bulk contribution does not depend on the choice of Ep and 
edge contributions for energies within A are finite and were calculated in 
[7]. Thus it suffices to consider Ep = —m. 

The required vacuum expectation values are 




sgn(m), if 7717 > 0, 
0, otherwise. 




Ige' 
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where 

dl dk 



•^'bulk ~ / -^ifclbulk 



27r2 ' 



iedge = / - ^) ifc ledge ^, 



and 

Jifc(x)|bulk = ^^lfc(x) O-^ U/fc(x), 
Jt(x)|edge = f/l(x)a^'f/fc(x). 

Here the variable E in n^fc is the negative roots of k'^ + P + m? . 

Note that j;\|buik and ledge vanish identically. On the other hand eqs. 
(4, 5) yield 



(7) i4(x)|buik = ;^-^3?(^^e 

where 

g = m — E + -y{k — il), 
f = {k-il)g\ 

The first term on the r.h.s. of eq. (7) is an odd function of A;, so that 
its contribution to the vacuum expectation value of vanishes for any 
symmetric regularization. For the second term we also will find a partial 
cancelation. 

We substitute k G (-C)0,oo) hy v = (^■■^'^s%nh(k/a) ^ (0,00), which yields 
for negative E 



(8) 

Moreover 
Here 



dk dv 



gE 



dv. 



Dl = |(7 + - V3){v - V4), 



where 



BOUNDARY STATES AND EDGE CURRENTS FOR FREE FERMIONS 



and 

il -\- m 7 — 1 

\/P + m? 7 + 1 ' 
—il + m 



The two denominators Di and D2 are formally related by a reflection across 
the boundary, with (m,, 7) 1— ?• (— m,7~"'^). 

One obtains a partial fraction expansion 

f flh 

--^ = {Pl{v) + P2{V) + P3(^) + Pa{v)] dv 

g E 

with 

a 



P2(^;) = --v 



a _2 



To regularize the momentum integrals we have to introduce a cutoff. In a 
solid state context the underlying physics selects a stationary Lorentz frame, 
in which a symmetric cutoff in k is natural. For this cutoff the dv integrals 
of these two terms cancel against each other due to the symmetry k —k, 
V I— 7- 1/v. The remaining terms are 

Ps{v) =ill±lv~\ 

7-1 

47 _-! 

P^(^v) = -il ^(^_^;3) . 

7^ — 1 

For a symmetric cutoff A of the t;-integration one finds 

y«(P3(„).--)^ = -^l±i.'(.),nA. 



To calculate the P4 integral we need 



il + m 7 — 1 
\//2"Tm2 7 + 1 



dv = ln(A) - In 



1 + 7 m + il 



1 — 7 V m — il 



where the principal value of the logarithm has to be taken, so that 

1 + 7 



In 



1 + 7 m + il 



1 — 7 V m — il 



In 



m + il 
m — il 



+ In 



1-7 



Tii 0(7 
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To calculate the dl integral over the logarithmic part we use 




m + il\ , , ^ „ f°° , / m + il\ r,-, „ 
—— cos{2lx) dl= \ J , e^*^^ dl. 

m — il I J~oo \ V m — il I 



The argument of In has branching points at I = zizim. We may move the 
integration path towards the cut with branching point I = im. Since the 
discontinuity of the logarithm is in, this yields an elementary integral. Al- 
together we obtain 



(9) j'(x)buik 

'1 + 7 



1-7 



6'ix) 



2nij^ - 1) \2x^ X J 7r(72 - 1) 2x2 

For the contribution from the edge functions, we have according to (3) 

7^-1 
dk = dX. 

7"^ + 1 

We have to integrate over those A for which A > and E < —m. This yields 



■2( ^ _ 7 



^(72 - 1) 



Note that j^dge vanishes for 7 E (—1,0) and is non-singular for 7 > 1. 

For 72 > 1 both the bulk and the edge current densities decay alge- 
braically at large distance from the boundary. In their sum these terms 
cancel, an interesting effect which reflects the tangency of the edge state 
and bulk state manifolds in [k, E) space, which was discussed in the previ- 
ous section. 

In all cases the singular part of the vacuum expectation value has the 
form 

.2.^\ _ 1 7' + l.,,^^ . l7l 1 



(j2(x)\ . , = 1-^s'ix) In A 



'singular 27r 72 - 1 ^ ' 47r{-f'^ - 1) ^2 ' 

Remarkably, this equation is invariant under y reflection together with 7 1— t- 
—7"^, such that it remains true for m < and is in fact independent of m. 
The regular part of the expectation value is continuous in m and vanishes for 
m = 0. For infinite m the contribution becomes localized at the boundary 
and is given by the 5' terms in (9). It may be interesting to rederive this 
result from a Chern-Simons Lagrangian. 
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5 Conclusion 



A single free fermion yields a divergent edge current on a half plane. For 
N fermions with boundary conditions 7„ the cut-off dependent edge current 
vanishes for 

^ 2,1 



=1 

One still has a dipolar edge current 

N 



0. 



E 



7n 



1 



■In 



1+7,1 



1 - 7n. 



which in general does not vanish. In addition there is a l/x^ singularity, 
which yields an unphysical divergent particle transport in the neighborhood 
of the edge, unless 

TV 



E 

n=l 



I7n| 



1 



0. 



Since the bulk system and the boundary at x = are Lorentz invariant 
under boosts in the y direction, one should be able to act with Lorentz 
transformations on the boundary conditions. In the context of relativistc 
quantum mechanics this is straightforward. The action on the parameters 
7n can be determined from eq. (3) or more easily from eqs. (8, 9), according 
to which boundary states are characterized by a Lorentz invariant signature 
rjn and their fixed travel velocity 



27r^ 

l+7n 

This translates into a rapidity On given by tanh0„ = Vn, or equivalently by 

l+7n 



rjn exp Or, 



1 - 7ri 



to which the rapidity of the Lorentz boost is added. When one moves from 
one-particle states to field theory, control of divergencies requires a cutoff 
which breaks Lorentz invariance. It may or may not be restored, when this 
cutoff is taken to infinity. One would expect that a restoration of Loretz 
invariance implies CPT invariance. Due to the unbroken PT invariance this 
should yield invariance under charge conjugation and consequently an ab- 
sence of net boundary currents. In our case we will find an unexpectedly 
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subtlety, however. In terms of edge velocities and signatures the two diver- 
gency cancellation conditions can be rewritten as 

N N 

^ r/„ exp(e„6'„) = = ^ r/„, exp(-e„,6'„), 

n=l n=l 

where 

e„ = sgn't;^. 

For = 1 there is no solution, in agreement with the restriction to an even 
N coming from the bulk physics. For N = 2 one needs 72 = —Ji^, thus 
invariance under charge conjugation. More generally, pairs of edge states 
related by charge conjegation do not contribute to the breaking of Lorentz 
invariance or to net boundary currents, as expected. There is an unexpected 
remainder of Lorentz invariance in more general situation, however. When 
all Vn of the remaining states have the same sign, the cancellation conditions 
are invariant under sufficiently small boosts which do not change the sign 
of edge velocities. When f„ of different signs occurs, no aspect of Lorentz 
invariance is recovered when the cutoff is removed. In any case, breaking of 
Lorentz invariance by edge states does not seem to be an argument against 
the realization of such systems in solid state physics, where this invariance 
is only a low energy phenomenon. 

It will be interesting to see to what extent the 7 values can be controled 
experimentally. Systems without time reversal invariance would be partic- 
ularly interesting, but are certainly very difficult to realize. 
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